/ Celestin Apprentice 4 / Apprentice-Release4.iso / Languages / RLaB 1.18c / manual / eigs.tex

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 5c 73 65 63 74 69 6f 6e | 7b 65 69 67 73 20 2d 2d |\section|{eigs --|
|00000010| 2d 20 53 79 6d 6d 65 74 | 72 69 63 20 45 69 67 65 |- Symmet|ric Eige|
|00000020| 6e 20 44 65 63 6f 6d 70 | 6f 73 69 74 69 6f 6e 7d |n Decomp|osition}|
|00000030| 0d 5c 6c 61 62 65 6c 7b | 66 6e 3a 65 69 67 73 7d |.\label{|fn:eigs}|
|00000040| 0d 5c 69 6e 64 65 78 7b | 65 69 67 73 28 29 7d 0d |.\index{|eigs()}.|
|00000050| 5c 69 6e 64 65 78 7b 65 | 69 67 65 6e 20 64 65 63 |\index{e|igen dec|
|00000060| 6f 6d 70 6f 73 69 74 69 | 6f 6e 7d 0d 5c 69 6e 64 |ompositi|on}.\ind|
|00000070| 65 78 7b 73 79 6d 6d 65 | 74 72 69 63 20 65 69 67 |ex{symme|tric eig|
|00000080| 65 6e 20 64 65 63 6f 6d | 70 6f 73 69 74 69 6f 6e |en decom|position|
|00000090| 7d 0d 5c 62 65 67 69 6e | 7b 72 61 69 6c 7d 0d 45 |}.\begin|{rail}.E|
|000000a0| 69 67 65 6e 53 20 3a 20 | 27 65 69 67 73 28 27 20 |igenS : |'eigs(' |
|000000b0| 28 20 4d 61 74 72 69 78 | 20 7c 20 4d 61 74 72 69 |( Matrix| | Matri|
|000000c0| 78 20 27 2c 27 20 4d 61 | 74 72 69 78 20 29 20 27 |x ',' Ma|trix ) '|
|000000d0| 29 27 20 3b 0d 5c 65 6e | 64 7b 72 61 69 6c 7d 0d |)' ;.\en|d{rail}.|
|000000e0| 0d 5c 76 65 72 62 2b 65 | 69 67 73 2b 20 77 69 74 |.\verb+e|igs+ wit|
|000000f0| 68 20 61 20 73 69 6e 67 | 6c 65 20 73 71 75 61 72 |h a sing|le squar|
|00000100| 65 20 61 72 67 75 6d 65 | 6e 74 20 6d 61 74 72 69 |e argume|nt matri|
|00000110| 78 20 63 6f 6d 70 75 74 | 65 73 20 74 68 61 74 20 |x comput|es that |
|00000120| 6d 61 74 72 69 78 27 73 | 0d 65 69 67 65 6e 76 65 |matrix's|.eigenve|
|00000130| 63 74 6f 72 73 20 61 6e | 64 20 65 69 67 65 6e 76 |ctors an|d eigenv|
|00000140| 61 6c 75 65 73 2e 20 54 | 68 65 20 72 65 73 75 6c |alues. T|he resul|
|00000150| 74 73 20 61 72 65 20 72 | 65 74 75 72 6e 65 64 20 |ts are r|eturned |
|00000160| 61 73 20 61 20 6c 69 73 | 74 20 77 69 74 68 0d 65 |as a lis|t with.e|
|00000170| 6c 65 6d 65 6e 74 73 20 | 5c 76 65 72 62 2b 76 61 |lements |\verb+va|
|00000180| 6c 2b 20 61 6e 64 20 5c | 76 65 72 62 2b 76 65 63 |l+ and \|verb+vec|
|00000190| 2b 20 77 68 69 63 68 20 | 61 72 65 20 74 68 65 20 |+ which |are the |
|000001a0| 65 69 67 65 6e 76 61 6c | 75 65 73 20 61 6e 64 0d |eigenval|ues and.|
|000001b0| 72 69 67 68 74 20 65 69 | 67 65 6e 76 65 63 74 6f |right ei|genvecto|
|000001c0| 72 73 20 72 65 73 70 65 | 63 74 69 76 65 6c 79 2e |rs respe|ctively.|
|000001d0| 20 49 6e 20 67 65 6e 65 | 72 61 6c 2c 20 74 68 65 | In gene|ral, the|
|000001e0| 73 65 20 77 69 6c 6c 20 | 62 65 0d 63 6f 6d 70 6c |se will |be.compl|
|000001f0| 65 78 20 71 75 61 6e 74 | 69 74 69 65 73 2e 20 20 |ex quant|ities.  |
|00000200| 0d 0d 5c 76 65 72 62 2b | 65 69 67 73 2b 20 77 69 |..\verb+|eigs+ wi|
|00000210| 74 68 20 74 77 6f 20 61 | 72 67 75 6d 65 6e 74 73 |th two a|rguments|
|00000220| 20 63 6f 6d 70 75 74 65 | 73 20 74 68 65 20 65 69 | compute|s the ei|
|00000230| 67 65 6e 76 65 63 74 6f | 72 73 20 61 6e 64 20 0d |genvecto|rs and .|
|00000240| 65 69 67 65 6e 76 61 6c | 75 65 73 20 6f 66 20 74 |eigenval|ues of t|
|00000250| 68 65 20 72 65 61 6c 20 | 67 65 6e 65 72 61 6c 69 |he real |generali|
|00000260| 7a 65 64 20 73 79 6d 6d | 65 74 72 69 63 20 28 20 |zed symm|etric ( |
|00000270| 6f 72 20 63 6f 6d 70 6c | 65 78 20 67 65 6e 65 72 |or compl|ex gener|
|00000280| 61 6c 69 7a 65 64 0d 48 | 65 72 6d 69 74 69 61 6e |alized.H|ermitian|
|00000290| 29 20 64 65 66 69 6e 69 | 74 65 20 65 69 67 65 6e |) defini|te eigen|
|000002a0| 70 72 6f 62 6c 65 6d 2e | 20 54 68 65 20 65 69 67 |problem.| The eig|
|000002b0| 65 6e 76 61 6c 75 65 73 | 20 61 6e 64 20 65 69 67 |envalues| and eig|
|000002c0| 65 6e 76 65 63 74 6f 72 | 73 20 61 72 65 0d 72 65 |envector|s are.re|
|000002d0| 74 75 72 6e 65 64 20 69 | 6e 20 61 20 6c 69 73 74 |turned i|n a list|
|000002e0| 20 61 73 20 66 6f 72 20 | 74 68 65 20 66 69 72 73 | as for |the firs|
|000002f0| 74 20 66 6f 72 6d 2e 20 | 20 49 66 20 74 68 65 20 |t form. | If the |
|00000300| 74 77 6f 20 61 72 67 75 | 6d 65 6e 74 73 20 61 72 |two argu|ments ar|
|00000310| 65 0d 24 41 24 20 61 6e | 64 20 24 42 24 2c 20 74 |e.$A$ an|d $B$, t|
|00000320| 68 65 20 65 69 67 65 6e | 70 72 6f 62 6c 65 6d 20 |he eigen|problem |
|00000330| 69 73 20 64 65 66 69 6e | 65 64 20 62 79 0d 0d 5c |is defin|ed by..\|
|00000340| 5b 20 41 78 20 3d 20 5c | 6c 61 6d 62 64 61 20 42 |[ Ax = \|lambda B|
|00000350| 20 78 20 5c 5d 0d 0d 54 | 68 69 73 20 72 6f 75 74 | x \]..T|his rout|
|00000360| 69 6e 65 20 66 6f 72 63 | 65 73 20 75 73 65 20 6f |ine forc|es use o|
|00000370| 66 20 74 68 65 20 73 79 | 6d 6d 65 74 72 69 63 20 |f the sy|mmetric |
|00000380| 65 69 67 65 6e 73 6f 6c | 76 65 72 20 2d 2d 2d 20 |eigensol|ver --- |
|00000390| 6e 6f 20 63 68 65 63 6b | 69 6e 67 0d 69 73 20 64 |no check|ing.is d|
|000003a0| 6f 6e 65 2e 0d 0d 5c 70 | 61 72 61 67 72 61 70 68 |one...\p|aragraph|
|000003b0| 7b 45 78 61 6d 70 6c 65 | 7d 0d 5c 69 6e 64 65 78 |{Example|}.\index|
|000003c0| 7b 65 78 61 6d 70 6c 65 | 20 6f 66 20 65 69 67 73 |{example| of eigs|
|000003d0| 28 29 7d 0d 5c 69 6e 64 | 65 78 7b 65 69 67 73 28 |()}.\ind|ex{eigs(|
|000003e0| 29 21 65 78 61 6d 70 6c | 65 7d 0d 5c 62 65 67 69 |)!exampl|e}.\begi|
|000003f0| 6e 7b 76 65 72 62 61 74 | 69 6d 7d 0d 3e 20 42 20 |n{verbat|im}.> B |
|00000400| 3d 20 72 61 6e 64 28 33 | 2c 33 29 0d 20 42 20 3d |= rand(3|,3). B =|
|00000410| 0d 20 20 20 20 20 30 2e | 35 37 20 20 20 20 20 20 |.     0.|57      |
|00000420| 30 2e 38 33 31 20 20 20 | 20 20 20 30 2e 34 36 32 |0.831   |   0.462|
|00000430| 20 20 0d 20 20 20 20 30 | 2e 36 38 33 20 20 20 20 |  .    0|.683    |
|00000440| 20 30 2e 30 37 33 35 20 | 20 20 20 20 20 30 2e 34 | 0.0735 |     0.4|
|00000450| 33 35 20 20 0d 20 20 20 | 20 30 2e 33 39 36 20 20 |35  .   | 0.396  |
|00000460| 20 20 20 20 30 2e 34 36 | 34 20 20 20 20 20 20 20 |    0.46|4       |
|00000470| 30 2e 36 33 20 20 0d 3e | 20 43 20 3d 20 73 79 6d |0.63  .>| C = sym|
|00000480| 6d 28 42 29 0d 20 43 20 | 3d 0d 20 20 20 20 20 30 |m(B). C |=.     0|
|00000490| 2e 35 37 20 20 20 20 20 | 20 30 2e 37 35 37 20 20 |.57     | 0.757  |
|000004a0| 20 20 20 20 30 2e 34 32 | 39 20 20 0d 20 20 20 20 |    0.42|9  .    |
|000004b0| 30 2e 37 35 37 20 20 20 | 20 20 30 2e 30 37 33 35 |0.757   |  0.0735|
|000004c0| 20 20 20 20 20 20 20 30 | 2e 34 35 20 20 0d 20 20 |       0|.45  .  |
|000004d0| 20 20 30 2e 34 32 39 20 | 20 20 20 20 20 20 30 2e |  0.429 |      0.|
|000004e0| 34 35 20 20 20 20 20 20 | 20 30 2e 36 33 20 20 0d |45      | 0.63  .|
|000004f0| 3e 20 65 69 67 65 6e 53 | 20 3d 20 65 69 67 73 28 |> eigenS| = eigs(|
|00000500| 43 29 0d 20 65 69 67 65 | 6e 53 20 3d 0d 20 20 20 |C). eige|nS =.   |
|00000510| 76 61 6c 20 20 20 20 20 | 20 20 20 20 20 76 65 63 |val     |     vec|
|00000520| 20 20 20 20 20 20 20 20 | 20 20 0d 3e 20 43 2a 20 |        |  .> C* |
|00000530| 65 69 67 65 6e 53 2e 76 | 65 63 5b 3b 31 5d 20 2d |eigenS.v|ec[;1] -|
|00000540| 20 65 69 67 65 6e 53 2e | 76 61 6c 5b 31 5d 20 2a | eigenS.|val[1] *|
|00000550| 20 65 69 67 65 6e 53 2e | 76 65 63 5b 3b 31 5d 0d | eigenS.|vec[;1].|
|00000560| 20 32 2e 37 38 65 2d 31 | 36 20 20 0d 20 35 2e 35 | 2.78e-1|6  . 5.5|
|00000570| 35 65 2d 31 37 20 20 0d | 20 32 2e 32 32 65 2d 31 |5e-17  .| 2.22e-1|
|00000580| 36 20 20 0d 3e 20 43 2a | 20 65 69 67 65 6e 53 2e |6  .> C*| eigenS.|
|00000590| 76 65 63 5b 3b 32 5d 20 | 2d 20 65 69 67 65 6e 53 |vec[;2] |- eigenS|
|000005a0| 2e 76 61 6c 5b 32 5d 20 | 2a 20 65 69 67 65 6e 53 |.val[2] |* eigenS|
|000005b0| 2e 76 65 63 5b 3b 32 5d | 0d 2d 31 2e 33 39 65 2d |.vec[;2]|.-1.39e-|
|000005c0| 31 36 20 20 0d 2d 36 2e | 32 35 65 2d 31 37 20 20 |16  .-6.|25e-17  |
|000005d0| 0d 2d 31 2e 36 37 65 2d | 31 36 20 20 0d 3e 20 43 |.-1.67e-|16  .> C|
|000005e0| 2a 20 65 69 67 65 6e 53 | 2e 76 65 63 5b 3b 33 5d |* eigenS|.vec[;3]|
|000005f0| 20 2d 20 65 69 67 65 6e | 53 2e 76 61 6c 5b 33 5d | - eigen|S.val[3]|
|00000600| 20 2a 20 65 69 67 65 6e | 53 2e 76 65 63 5b 3b 33 | * eigen|S.vec[;3|
|00000610| 5d 0d 2d 34 2e 34 34 65 | 2d 31 36 20 20 0d 20 32 |].-4.44e|-16  . 2|
|00000620| 2e 32 32 65 2d 31 36 20 | 20 0d 20 20 20 20 20 20 |.22e-16 | .      |
|00000630| 20 20 30 20 20 0d 5c 65 | 6e 64 7b 76 65 72 62 61 |  0  .\e|nd{verba|
|00000640| 74 69 6d 7d 0d 0d 53 65 | 65 20 61 6c 73 6f 3a 20 |tim}..Se|e also: |
|00000650| 50 61 67 65 7e 5c 70 61 | 67 65 72 65 66 7b 66 6e |Page~\pa|geref{fn|
|00000660| 3a 65 69 67 7d 20 5c 76 | 65 72 62 2b 65 69 67 28 |:eig} \v|erb+eig(|
|00000670| 29 2b 2c 0d 50 61 67 65 | 7e 5c 70 61 67 65 72 65 |)+,.Page|~\pagere|
|00000680| 66 7b 66 6e 3a 65 69 67 | 6e 7d 20 5c 76 65 72 62 |f{fn:eig|n} \verb|
|00000690| 2b 65 69 67 6e 28 29 2b | 2e 0d 0d 5c 6e 65 77 70 |+eign()+|...\newp|
|000006a0| 61 67 65 0d 0d 0d 0d 0d |                         |age.....|        |
+--------+-------------------------+-------------------------+--------+--------+